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Abstract 

We obtain closed-form expressions, in terms of the Faulhaber 
numbers, for the weak-field expansion coefficients of the two-loop 
Euler-Heisenberg effective Lagrangians in a magnetic or electric 
field. This follows from the observation that the magnetic world- 
line Green's function has a natural expansion in terms of the 
Faulhaber numbers. 



1 Introduction: The Euler-Heisenberg Lagrangian 
at one and two loops 

The Euler-Heisenberg Lagrangian |2] describes the effect of a virtual elec- 
tron - positron pair on an external Maxwell field in the one loop and constant 
field approximation. Its standard proper time representation is 
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Here T is the proper-time of the loop fermion, m its mass, and a, h are the 
two Maxwell invariants, related to E, B by a 2 — b 2 = B 2 — E 2 , ab = E ■ B. 
The superscript " (1)" stands for one loop. A similar representation exists for 
scalar QED 0E|: 
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The Lagrangians (jl.lj) . (jl.2j) historically provided the first examples for the 
concept of an effective Lagrangian, and moreover the first nonperturbative 
result in quantum field theory. See |3] for a review of their many applications 
and generalizations. 

The proper time integrals in these formulas can be done exactly in terms 
of certain special functions |3] . Alternatively, one can expand the integrands 
as power series in the field invariants, using the Taylor expansions 
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Here the Bm are the Bernoulli numbers. The terms in this expansion involv- 
ing N = 2n powers of the field contain the information on the low energy 
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limits of the N photon scattering amplitudes. Thus in the low energy limit 

one can obtain these amplitudes in closed form [5]. 

For a purely magnetic field eqs. (JEU), (JUU), flO|), (d) yield 
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An analysis of the coefficients shows that these series are divergent but Borel 
summable [SHE]- Moreover, the summability relates to the fact that the 
effective Lagrangian is real in the magnetic case. The corresponding series 
for the purely electric case is obtained by the replacement B 2 — ► — E 2 , which 
turns the alternating series into a non-alternating one. The non-alternating 
series is not Borel summable, which is an indication of the well-known fact 
that the electric effective Lagrangian has an imaginary part. Although this 
imaginary part is nonperturbative in nature, it is possible to calculate it from 
the expansion coefficients by an analysis of their large n behaviour, combined 
with a Borel dispersion relation [El IE] ■ 

This paper addresses the possibility of analyzing such weak-field expan- 
sions at higher loop orders in QED. The two loop corrections to the La- 
grangians (jl.lj) . (jf .2j) have first been considered by Ritus, who obtained them 
in terms of integral representations both in spinor EI] and scalar QED 
[TT] . Other representations for the same Lagrangians were later given in 
[T2*t lT3*t ITU ITH] . however they all involve the same apparently intractable 
type of double integrals. As a consequence, at the two loop level presently 



only the first few coefficients of the power series expansions for £s pin / scal (-B) 

are known 8], and the imaginary parts Im£^ n / scal (E) are known explicitly 
only in the leading weak field limit jTHl IE] (although Lebedev and Ritus 
succeeded in establishing the general structure of the imaginary part [To]). 

This should be contrasted with the case of a self-dual (euclidean) field, 
corresponding to external photon lines of definite helicity, where the two-loop 
effective Lagrangians can be obtained in closed form ^TJEEIGS]- This result 
allows one to extend the aforementioned calculation of the low energy limit of 
the N photon amplitudes to the two loop level for one particular component 
of the photon S matrix, the so-called MHV amplitudes [THj . 

For general fields, it seems difficult to make further progress at the two 
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loop level without having a closed form expression for the weak field ex- 
pansion coefficients. It is the purpose of the present paper to derive such 
expressions for the purely magnetic (or purely electric) cases. Thus, our goal 
is to obtain the coefficients of the weak field expansion of the (renormalized) 
magnetic two-loop effective Euler-Heisenberg Lagrangians £s pi n/scai(-^)> 

- (s)>i>s>)(sr 2 > w 

n=0 

n=0 

Note the difference in normalization between the scalar and spinor QED 
cases, which takes into account the global factor of —2 from statistics and 
degrees of freedom. 

2 Two loop expansion coefficients for spinor 
QED 

We start with the following integral representation of the on-shell renormal- 
ized Lagrangian, which was obtained in [Tlj : 

4pin( 5 ) = 4pin,main( 5 ) + 4pin,<5m( 5 )- (2' 1 ) 

Here the "main part" is given by the following two parameter integral, 
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where z = eBT and 
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The integrand involves the so-called worldline Green's function G(u) and the 
"magnetic" Green's function G z (z,u), as well as their w-derivatives (denoted 
as G and G z , respectively): 

= u(l — u) , 
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The function L(z, u) is essentially the two-loop integrand before renormaliza- 
tion. The term Lq2{z,u) removes the order z°,z 2 terms, which implements 
the renormalization of the charge and the field, and the removal of the vac- 
uum energy. (In the following we will not always make this subtraction 
explicit.) The other subtraction term involving relates to mass renor- 
malization, which is necessary starting at the two loop level. In the worldline 
formalism the need for this mass renormalization subtraction can be recog- 
nized from the appearance of singularities at u = and u—1 ■ Those 
singular terms can be absorbed into a mass shift of the one-loop lagrangian, 
^ m °am^spin(-^); however this leaves a finite remainder C^ in Sm (B), 
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The weak field expansion coefficients of C^ in Sm (B) are again easy to obtain 
using the Taylor expansion (jl.3j) . One finds 
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However, the main integral term ()2.2|) is much more difficult to expand. 
A brute-force expansion was done in [Sj, but no closed-form expression was 
obtained. The technical challenge is to expand the integrand (excluding the 

m 2 

t~TB 2 factor) in a series in z, in such a way that the u integrals can also 
be done easily. This is complicated by the fact that the magnetic Green's 
function G z (z, u) couples these two parametric variables in a nontrivial man- 
ner. The main observation of this paper is that there exists an expansion of 
G z (z,u) which decouples the variables in an elegant manner. Remarkably, 
G z (z, u) is essentially the generating function of the so-called Faulhaber num- 
bers of number theory [20J. Namely, 
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Here we have defined 

f(m,k)= f(m,k) (2.8) 

(2m + 1)! 

with f(m, k) the Faulhaber numbers. These numbers can in turn be written 
as linear combinations of the Bernoulli numbers j2Hl : 
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Given such an expansion, the u integrals are now trivial as they involve 
powers of the free Green's function G(u) = u(l — u). Using the standard 
Euler beta function integral we define 

j^du[u{l-u)] n = {2 f +l)l =/?(n + l). (2.10) 
To illustrate the Faulhaber expansion (|2.7jl . consider the first few terms: 
G z = G + (l - z coth(z)) G 2 + 0(G 3 ) = G - l -G 2 z 2 + 0{z\ G 2 ) . (2.11) 

This suggests rearranging L(z,u) as a series in AG/G where AG = G z — G 
is the difference between the magnetic and free worldline Green's functions. 
After simple manipulations, this leads to 
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The u derivatives appearing in the last term can be removed by an integration 
by parts, yielding 
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Performing the subtraction, implied in ()2.2|) . of all terms O(^), we finally 
obtain 
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Here ()2.7|) was used to obtain the ^ subtraction for the last term. We now 
expand the trigonometric functions using (jl.Hj) and 
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This leads directly to the following closed-form expression for the expansion 
coefficients: 
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The full two-loop expansion coefficients in (|1.7|) are now given in closed-form 
by combining ()2.6|) and ()2.17|) : 
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3 Expansion coefficients for scalar QED 

The case of scalar QED can be treated completely analogously, starting from 
any of the various representations for cfJ^B) given in [T3J El UH] • We will 
give here only the final result for the weak field expansion coefficients (as 
defined by IfTTBjl ): 
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4 Conclusions 

We have found closed-form expressions for the weak field expansion coef- 
ficients of the two loop corrections to the renormalized Euler-Heisenberg 
Lagrangians in a purely magnetic (or purely electric) field. As a check, we 
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have verified that eqs. (|2.19|) and (|3.1|) indeed reproduce the known low order 
coefficients in these expansions |T3J |8| : 
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Although the closed-form formulae ()2.19)1 and (|3.1j) are significantly more 
complicated than the corresponding ones for the case of a self-dual field j!8( 
H§] . their structure is still similar insofar as they can, using (|2.9j) . be written 
in terms of folded sums of Bernoulli numbers with factorial coefficients. Of 
course, it is quite possible that these formulas can still be simplified. 

We note that the two-loop expansion coefficients are still rational num- 
bers, after extracting a factor of (a/n) 2 , just as the one-loop coefficients are 
rational after extracting a factor of (a/n). A question of obvious interest 
is whether this property persists to higher loop orders. Based on a com- 
parison with what is known about the coefficients of the QED (3 - functions 
[2H E2 E31 EH E5] we consider it likely that rationality will be found to hold 
at least for the quenched (order Nf) contributions to the Euler-Heisenberg 
Lagrangians at arbitrary loop order. 

Finally, since in the worldline formalism the magnetic Green's function 
G z is the basic ingredient appearing in the integral representations for all 
processes involving constant magnetic fields f3JEI] 5 we expect the Faulhaber 
expansion (J2.7)) to become useful for other calculations of this type, including 
possibly higher loop orders. 
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